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Abstract. - A new behavior type of unidirectionally coupled chaotic oscillators near the gen- 
eralized synchronization transition has been detected. It has been shown that the generalized 
synchronization appearance is preceded by the intermitted behavior: close to threshold param- 
eter value the coupled chaotic systems demonstrate the generalized synchronization most of the 
time, but there are time intervals during which the synchronized oscillations are interrupted by 
non-synchronous bursts. This type of the system behavior has been called intermitted general- 
ized synchronization (IGS) by analogy with intermitted lag synchronization (ILS) [Phys. Rev. 
E 62, 7497 (2000)]. 



Synchronization [1,2] of interacting chaotic oscillators is one of the fundamental phenomena 
of nonlinear dynamics. Recently several types of chaotic synchronization have been observed in 
coupled nonlinear oscillators. These are phase synchronization (PS) [3-5], lag synchronization 
(LS) [6-8], complete synchronization (CS) [9-12] and generalized synchronization (GS) [13- 
15]. It is important to note that GS may also takes place in the non-oscillatory chaotic systems 
(see [16] for detail). All synchronization types are interrelated (see for detail [7,13,17]), but 
the relationship between them is not completely clarified yet. In recent works [18-21] it has 
been shown that all these synchronization types may be considered from the common point 
of view as different manifestations of one universal phenomenon. 

It has also been found that onsets of phase and lag synchronization types are preceded by 
intermittent behavior. Close to the threshold parameter value the coupled chaotic systems 
demonstrate synchronized dynamics most of the time, but there are time intervals during 
which the synchronized oscillations are interrupted by non-synchronous behavior. These pre- 
transitional intermittencies have been described in [6,22,23] for the case of lag synchronization 
and in [24-26] for phase synchronization, respectively. 

Due to the existence of unifying framework of coupled chaotic oscillators synchronization 
one can also expect the intermitted behavior at the threshold of the generalized synchroniza- 
tion appearance. In this work we consider the behavior of coupled chaotic oscillators close 
to the coupling parameter value corresponding to the onset of generalized synchronization 
regime. As it will be shown below the generalized synchronization appearance is also pre- 
ceded by the intermitted behavior in the same way as the phase synchronization and lag 
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Fig. 1 - The dependence of difference between the coordinates of response and auxiliary systems 
Mi (t) — Vi (t) on time. The coupling parameter between drive and response systems has been selected 
as e = 0.106 



synchronization regimes do. This type of system behavior we called intermitted generalized 
synchronization (IGS) by analogy with intermitted lag synchronization (see [22]). 

The generalized synchronization (GS) [13-15] regime introduced for drive x(t) and response 
u(t) systems means that there is some functional relation between unidirectionally coupled 
chaotic oscillators, i.e. u(t) — F[x(i)]. This functional relation F[-] can be very complicated, 
but there are several methods to detect the synchronized behavior of coupled chaotic oscilla- 
tors, e.g., the method of nearest neighbors [13,27]. In our work we have used the auxiliary 
system approach proposed in [28] . In this case we consider the dynamics of drive 

*(*) = H(x(t)) (1) 

and response 

u(t) = G(u(t),g,x(t)) (2) 

systems. The vector g characterizes the coupling of the drive x(t) and response u(t) systems. 
At the same time we also consider the dynamics of the auxiliary system 

v(t) = G(v(t),g,x(i)) (3) 

which is identical to the response system u(t) but starts with the other initial conditions, 
i.e., u(io) 7^ v (^o)- In the absence of generalized synchronization between the drive x(i) and 
response u(t) systems, the phase trajectories of the response u(i) and auxiliary v(t) systems 
share the same chaotic attractor but are otherwise unrelated. In the case of generalized 
synchronization the orbits of response u(t) and auxiliary v(t) systems become identical after 
the transient dies out due to the generalized synchronization relations u(t) — F[x(i)] and 
v(t) = F[x(i)]. It is obvious that in the case of generalized synchronization the condition 
u(i) = v(i) should be satisfied and the identity of response and auxiliary systems is a much 
simpler criterion to test the presence of generalized synchronization rather than detection of 
unknown functional relationship F[-]. 

For the coupling parameter values close to the onset of generalized synchronization the 
drive and response systems demonstrate the regime of intermitted generalized synchronization. 
In this case the generalized synchronization regime between drive and response systems has 
been observed most of the time and at that time the vector states of response u(t) and 
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Fig. 2 - The statistical distribution of laminar phases and its approximation N(t) ~ t~ 3//2 in log-log 
scale. The coupling parameter between drive and response systems has been selected as e — 0.106, 
the value of the threshold A = 0.1 

auxiliary v(£) systems coincide with each other. Nevertheless, there are time intervals when 
the drive and response systems behave unrelatedly. Therefore, during these time intervals 
the vector states of response u(t) and auxiliary v(t) systems differ from each other. If we 
plot the difference between vector states of response and auxiliary systems |u(t) — v(t)| versus 
time variable t than the graph will be obtained where the laminar phases corresponding to 
the generalized synchronous behavior alternate with the turbulent ones when the generalized 
synchronization is not detected. In other words, in the case of IGS for some time intervals 
one can detect the presence of functional relation F[-] between the drive and response systems 
whereas for the other time intervals such functional relation does not exist. 

Let us consider the dynamics of two unidirectionally coupled Rossler systems 



where x = (xi,X2,X3) T and u = (uj., ua,U3) T are the vector states of drive and response 
systems, respectively. The control parameter values have been chosen as u>i — 0.99, u)<z — 0.95, 
a = 0.15, p = 0.2, c = 10, the parameter e characterizes the coupling strength between 
the considered systems. At the same time the behavior of auxiliary system v(t) described 
also by equation |@J should be considered in order to detect the presence of the generalized 
synchronization regime. As it follows from our investigation PS arises at sps ~ 0.116, GS 
takes place at e at e c w 0.11, and LS can be observed for s > Els ~ 0.36. So, GS for the 
coupling strength exceeding slightly the critical value e c is the week synchronization [15, 16] 
and, therefore, the relationship F[-] demonstrates fractal properties. For e > e^s GS becomes 
the strong synchronization as well as F[-] begins being smooth. In this case the behavior of 
the drive and response systems becomes identical but shifted in time versus each other. Close 
to the onset of the GS the regime of IGS has been detected. The difference between the 
coordinates of response and auxiliary Rossler systems has been shown in Fig. ^ 

The time intervals where the difference U\{t) — Vi(t) is close to zero correspond to the 
laminar phases of coupled oscillator dynamics when the presence of functional relation F[-] 
between drive and response systems can be detected. Alternatively, the irregular bursts in- 



Xl = -LUlX 2 - x 3 
X2 = UJ1X1 + aX2 
X 3 =p + X 3 (Xi - 



and 



ill = -u 2 u 2 - u 3 + e{x\ - tti) 
U2 = uj 2 ui + au 2 
u 3 =p + u 3 (m - c), 



(4) 



4 



EUROPHYSICS LETTERS 



200 




100 I . . . . I 

0.001 0.002 0.005 s-e c 

Fig. 3 - The log-log plot of the mean laminar phase length (r) vs. the deviation (e c — e) of coupling 
strength from the critical value and its approximation by (r) ~ (e c — e) _1 . The power low with 
critical exponent —1 has been observed 

dicate breaking of the relationship F[-]. It is interesting to note that amplitude of irregular 
bursts is large enough and comparable with the size of chaotic attractor of Rossler system. 

This type of behavior demonstrates the features of on-off intermittency [29, 30] which 
usually appears at the desynchronization of two coupled identical chaotic oscillators. To 
analyze the statistical properties of IGS we numerically compute the distribution N(r) of 
the laminar phases length r. We suppose that the current state is the laminar phase if the 
difference between variables U\(t) and V\(t) is less than A = 0.1, otherwise the current state 
assumes to be an irregular burst. 

Distribution of the laminar phases length has been shown in Fig. [21 One can see that this 
distribution is close to the power law with exponent n = —3/2. It should be noticed that this 
result does not sensitively depend on the value of threshold A. 

The other criteria of the on-off intermittency presence is the dependence of the mean 
laminar phase length (r) on the deviation of coupling strength e from the critical value e c . 
Fig. shows the numerically determined mean laminar phase length (r) versus deviation 
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Fig. 4 - Circuit realization of two unidirectionally coupled Chua's oscillators 
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Fig. 5 - The dependence of difference ui (t) — t>i (t) between the coordinates of response and auxiliary 
Chua's circuits on dimensionless time £. The coupling parameter between drive and response systems 
has been selected as e = 0.03105, the control parameters are a x = 2.78, a u = 2.89, 7 = 3.00 

(e c — e). From this figure one can see the universal power low (r) ~ (e c — e) with critical 
exponent k = — 1. Therefore, the IGS appears to be an on-off intermittency. It is important 
to note, that the ILS is also the on-off intermittency regime (see for detail [23]). 

It is important to note that the attractor of the considered Rossler system is coherent 
for the selected values of the control parameters. At the same time we have observed the 
IGS presence for the different values of the control parameters when the system behavior 
is characterized by the funnel attractor (see, e.g., parameter values pointed in [5]). The 
character of chaotic attractor does not seem to be important for IGS as well as the parameter 
mismatch influence: IGS may be observed both for small and large mistuning of the oscillators 
parameters. 

The intermitted generalized synchronization regime seems to appear in the different cou- 
pled dynamical systems. To examine this assumption we have considered the dynamics of 
two unidirectionally coupled Chua's circuits [31,32] with piecewise linear characteristic (see 
Fig-EJ- The typical property of the considered Chua's circuit is the presence of two distinct 
characteristic time scales. Therefore, this system may demonstrate periodical, quasi-periodical 
and chaotic oscillations depending on the values of the control parameters. The behavior of 
the considered Chua's circuits in the autonomous regime is chaotic due to the selection of the 
control parameter values. 

The behavior of the drive Chua's circuit is described by 




X\) + x 3 ) , 



(5) 



and the equations of the response system are the following 



~f(u2 - 

1 £ 
(f(u-2 - Ul) + U3) + -[X2 — U 2 ), 

7 7 
7U2. 



(6) 



6 



EUROPHYSICS LETTERS 



The variables {x,u}i = V^ u / 'E and {x,u}2 = V c ^ /E are dimensionless voltages on ca- 
pacitors Ci' u and C2 of the first (drive) and the second (response) oscillators, respectively. 
The variables {x, u}^ = i x £ u / 'I are the dimensionless currents through the inductances L in 
the drive and response circuits (see Fig.0]for detail). The parameters E and / are the normal- 
ization factors. Dimensionless control parameters are a x ,u — C-2,jC\ ,u and 7 = -^-yC^/L; 
£ = tj \JLCi is dimensionless time. The function 

/W = -^+|(^ ± ^)(|. + H-k-l|) ! (7) 
mi 2 \ mi J 

is the dimensionless piecewise linear voltage-current characteristic of nonlinear element N, 
where uiq and mi are the conductivities of the corresponding branches of voltage-current 
characteristic. The ratio mo/mi has been selected as 1/2. The coupling parameter s = 1/ ' Rmi 
determines the influence of the drive Chua's circuit on the response one. 

In Fig. [5] the dependence of difference between the dimensionless voltages on capacitors 
Ci' v of response and auxiliary system on dimensionless time £ has been shown. One can 
see the presence of intermitted generalized synchronization between the drive and response 
Chua's circuits. 

Let us briefly discuss the nature of IGS. For this purpose we consider the extended phase 
space Dffilffi A consisting of the phase space of the drive (IS), response (R) and auxiliary A 
systems. We can also consider the subspaces D © R and 1® A of O © R © A. The presence 
of GS between the drive and response systems implies the complete synchronization between 
the response and auxiliary systems. As the criterion of GS one can consider the negativity of 
the values of the conditional Lyapunov exponents (CLEs) [10,15] in the OffiR phase subspace 
while the criterion of CS is the negativity of the transversal Lyapunov exponents (TLEs) [2] 
in the R © A subspace. Following the ideas of [15] one can draw a conclusion that CLEs (in 
D © R) and TLEs (in R © A) coincide with each other. Therefore, when GS is destroyed in 
B©R with coupling parameter decreasing, the blowout bifurcation [33,34] takes place in the 
M © A subspace. At the same time it is well known that on-off intermittency takes place 
below the coupling parameter value corresponding to the blowout bifurcation. So, the on-off 
intermittency has to appear in the ED © R subspace below the threshold of the GS arising. 

Taking into account this idea one can draw another conclusion concerning the noise in- 
fluence. As the noise results in the bubbling in the case of the identical chaotic oscillators 
for the coupling parameter values above the point of blowout bifurcation, the similar effects 
will be probably observed for GS. It means that the intermitted behavior may be observed in 
the experiments for the coupling parameter values corresponding to the GS regime due to the 
noise influence. 

In conclusion, we have detected the new behavior type of unidirectionally coupled chaotic 
oscillators near to the onset of generalized synchronization. This type of behavior has been 
called intermitted generalized synchronization. It should be noted that the similar effect has 
been observed for the unidirectionally coupled time-delayed chaotic oscillators [35]. As GS 
phenomenon could also be found in the non-oscillatory chaotic systems (see, for detail the 
monographic review [16] dealing with synchronization phenomena in chaotic systems) IGS 
can also be found in such systems below the threshold of GS arising, correspondingly. 

The intermitted generalized synchronization seems to be a common phenomenon in nonlin- 
ear systems, and two unidirectionally coupled chaotic oscillators usually undergo the transition 
to the generalized synchronization through the on-off intermittency. The intermitted behavior 
of coupled chaotic oscillators close to onset of PS, LS and GS is likely to be the manifestation 
of universalities underlying the synchronization phenomenon. 
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